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Abstract. We discuss how to use the Rees-Sciama (RS) effect associated with merging clusters of galaxies to measure their 
kinematic properties. In a previous work (Rubiiio-Martin et al. 2004), the morphology and symmetries of the effect were 
examined by means of a simplified model. Here, we use realistic N-body simulations to better describe the effect, and to 
confirm that the signal has a characteristic quadrupole structure. From the amplitude of the signal obtained, we conclude that it 
is necessary to combine several cluster mergers in order to achieve a detection. Using the extended Press-Schechter formalism, 
we characterized the expected distribution of the parameters describing the mergers, and we used these results to generate 
realistic mock catalogues of cluster mergers. To optimize t he ext raction of the RS signal, we developed an extension of the 
spatial filtering method described in Haehnelt & Teemark ('1996"). This extended filter has a general definition, so it can be 
applied in many other fields, such as gravitational tensing of the CMB or lensing of background galaxies. It has been applied 
to our mock catalogues, and we show that with the announced sensitivities of future experiments like the Atacama Cosmology 
Telescope (ACT), the South Pole Telescope (SPT) or the Atacama Large Millimeter Array (ALMA), a detection of the signal 
will be possible if we consider of the order of 1, 000 cluster mergers. 



1. Introduction 

Observations of the fluctuations of the cosmic microwave back- 
ground (CMB) can give us information about the formation 
of structures in the Universe, because the evolution of the 
gravitational potentials leaves its imprint on the CMB pho- 
tons as they travel along them. This physical effect is usu- 
ally split int o two terms. One is th e integrated Sach s-Wolfe 
effect (ISW; LSachs & Wolfl 1 19671: IHu & Sugivamal [l994) . 
which is produced by the linear evol ution of the potentials . 
The other is th e Rees-Sciama effect (RS; Rees & Sciama, 1968; 
iMartinez-Gonzalez et al.. J390; Seliak. .1996 ), in which the 
density contrast producing the gravitational potential is in its 
non- linear regime. 

In this work we discuss the RS effect associated with 
the non-li near regime of ga laxy clusters mergers. In a previ- 
ous work jRubino-Martm et al. (2004), hereafter Paper I), this 
problem was examined using a simplified model of the physics 
of the merger. The amplitude, morphology and symmetries of 
the effect were discussed and characterized in terms of the 
physical parameters of the merger. However, the RS effect is 
not the only signal being generated by clusters. The presence 
of hot gas in the Intra-Cluster Medium (ICM) and the peculiar 
velocity of the cluster with respect to the Hubble flow introduce 
temperature fluctuations which are proportional to the inte- 



gral of pressure (thermal Sunyaev-Zel'dovich effect, [tSZ]) and 
the momentum density (kinematic Sunyaev-ZeFdovich effect, 
[kSZ]) along the Une of sig ht, (.Sunvaev & Zeldovicht 119721; 
ISunvaev & Zel'dovichlll98(]l) . In addition, galaxy clusters may 
host radio-galaxies or infrared galaxies whose flux might be 
relevant to the small-amplitude effects we are studying here. 
Further, the intrinsic CMB field and its deflection caused by the 
matter distribution in the cluster produce different anisotropics 
which must be accounted for. In Paper I, it was proposed to co- 
add coherently the RS signals from a sample of cluster mergers 
in order to achieve a statistical detection of the effect against 
these "foreground" anisotrpy signals. 

Here, we extend this work in three ways. First, we use hy- 
drodynamical simulations in order to compute realistic maps 
of the RS effect, and to show that the morphology of the signal 
has the characteristic quadrupole structure described in Paper 
I. We also show that N-body simulations contain an additional 
signal which originates during the collapse of the environmen- 
tal dark-matter in the potential well of the system, although the 
morphology of this signal is completely different from the one 
associated with the merging induced one. Second, we perform 
a detailed study of the expected distribution of the parameters 
describing a merger of galaxy clusters (masses of the compo- 
nents, distances and velocities), and we present semianalyti- 
cal expressions for the probability of finding a cluster merger 
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with a given parameter set. These equations can be used to pro- 
duce a realistic "catalogue" of cluster mergers, which we use 
in order to forecast the detectability of the effect. Finally, we 
investigated the prospects for the RS signal extraction from 
cluster mergers by future CMB experiments, like the South 
Pole Telescope' or the Atacama Cosmology Telescope^. To this 
end, we extended the filtering method described in Haehnelt & 
Tegmark (1996). Our new filter has a general mathematical def- 
inition, so it could be applied to extract the signal in many other 
fields, e.g. gravitati onal lensing of the CMB ( Seljak^ .1996 : 
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2005a| background galaxies (e.g. 

2005b). We use this extended filter to separate 
the RS signal from the following three components: the primor- 
dial CMB, the gravitational lensing imprint, and the kinetic SZ 
effect. We show that with upcoming instruments, the RS sig- 
nal could be detected if we consider of the order of 1,000 ob- 
served cluster mergers. As we will see below, the main limiting 
"noise" component in such a measurement will be the lensed 
CMB. 



2. The Rees-Sciama effect 

In Paper I, the Rees-Sciama effect associated with mergers of 
galaxy clusters was discussed using an analytical model for 
the merger In addition, an analytical prescription to prepare 
Rees-Sciama maps from hydrodynamical simulations was in- 
troduced. Here, we implement this prescription and we present 
some RS maps resulting from numerical simulations. The mor- 
phology of the effect can then be directly compared to that pre- 
dicted in Paper I. 

In order to prepare a map of the RS effect from a hydrody- 
n amic simulation, we u s e the analytic recipe proposed in Paper 
I jRubino-Martm et al.1 1200 4) to compute the CMB tempera- 
ture change along the line-of-sight (LOS) direction n. 



(1) 



where g(r) is the total mass density and v(r) the mean velocity 
of matter. In this equation, the vectors are split into line-of- 
sight (LOS) parallel and LOS-perpendicular components, so 
r - irj_,z)- The axis along the line-of-sight is specified by 
h, with the positive values pointing towards the observer To 
be valid, this expression requires the gravitational fields to be 
weak (which is the case for galaxy clusters), and the matter 
distribution to be nearly stationary (i.e. v/c <K 1), so the mass 
density does not show an explicit dependence on time. The al- 
gorithm was integrated into the SMAC code^, which aims at 
deriving maps for the main galaxy cluster observables from nu- 
merical simulations. 

Due to its gravitational nature, the RS effect can also be re- 
lated to the deflection angle, yielding an expression very useful 
for analytical computations: 



-—(r^) = — siny8H„i ■ a{n) . 
T c 



(2) 



' SPT, see http://spt.uchicago.edu/| 

^ ACT, see http://www.hep.upenn.edu/~angelica/act/act.htmr 
^ (Simulated MAps Creator: http://dipastro.pd.as tro.it~cosmoj 



Here a is the gravitational lensing deflection field of the cluster, 
H,,x is the unitary transverse velocity vector and /? is the angle 
described by the cluster velocity v a nd the plane of the sky, as 
shown by Birkinshaw & Gulll ( Il983h . 

The dot product m,,± • «(«) has a dipolar structure elongated 
along and, according to Equation the RS anisotropy 
is proportional to the halo mass and to its transverse velocity. 
Due to these dependencies, forming galaxy clusters with two 
major merging haloes seem to be the best candidates for RS 
detection. First, because of the large masses and velocities in- 
volved, of the order of IO'^'Mq and l,000km/s respectively. 
Second, because the contributions of the merging components 
add constructively due to the converging cluster motions. 

A simple description of this scenario can be built by 
considering two dark matter haloes, having a NFW profile 
jNavarroeta'llll995h . in free fall. It provides a simple deriva- 
tion for the lensing and the RS effects, from which it is straight- 
forward to obtain the templates necessary for the data reduc- 
tion (Section|4j. The RS effect can be computed according to 
Equation ^ w here the defle ction angle can be computed ana- 
lytically ( Bartelman n. 1996). As it will be shown, the fact that 
the deflection field is axially symmetric on large scales justifies 
the assumption of the halos being spherical. 

As an example, we present in Figure Q the RS effect of a 
major merger event. In the left panel, the RS signal was com- 
puted with our simple analytical recipe, whereas in the right 
panel we used the numerica lly simulated cluster g24-i-200 de- 
scribed in iRasia et al ] ll2004 . This cluster consists of two merg- 
ing haloes, of masses Ml = 6xl0''*Mo//zandM2 = IO'^Mq/Zj, 
with an infall velocity of 1800A:ot/s, provided by G. Tormen. 
The field of view is 2 degrees on a side. This figure demon- 
strates that our simple analytical model can describe the RS 
signal generated in major cluster mergers with reasonable accu- 
racy. The main differences between our description and the out- 
put obtained from the numerical simulations are related to the 
presence of minor merging substructures within the main ob- 
ject, and the roughly spherical collapse of environmental dark 
matter (DM) onto the potential. 

The effect of minor mergers can be neglected, since their 
amplitude is small and their contribution vanishes after apply- 
ing our filter (see Section l4!2l below). The roughly spherical 
collapse of environmental DM has the shape of a negative cir- 
cular symmetric peak centered on the potential minimum, and 
adds coherently to the merger-induced RS signal. Also this ex- 
tra RS component is related to the dynamical status of DM and 
could be considered part of the signal we are searching for. 
The inclusion of this contribution in our procedure can be eas- 
ily implemented by simply changing the filter RS template (see 
Section In this case we would measure the full DM mo- 
mentum instead of only the DM momentum of the main merg- 
ers. Since we ignore this signal, our estimate of the detectability 
of the RS signal is conservative. 

3. Distribution of physical parameters describing 
the merger 

We investigate the probability distribution of the physical 
parameters describing a merger. We use the approach of 
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Fig. 1. A comparison between our analytical description of 
the RS signal generated in a major merger (left panel) and 
th e numerica l simulation g24+200 described in the text and 
in lRasia et alJ 1,2004.) (right panel). The field of view is 2 de- 
grees on a side. The two clusters are placed at redshift 0.3, 
have masses of Ml = 6 X 10''*Mo//i and M2 = IO'^Mq/Zz 
respectively, and they are falling onto each other with a rel- 
ative velocity of 1800 km/s. The analytical model represents 
the same merger of the numerical simulation, showing that our 
simple recipe provides a good description of the cluster merger 
The total contribution of minor mergers is negligible, as ex- 
plained in Section J4.2l i. The main difference is that the simu- 
lation includes the roughly spherical collapse of environmental 
DM onto the potential well. This collapse-induced signal has 
the shape of a negative peak with circular symmetry centered 
on the potential minimum. This contribution adds coherently to 
the signal of the two main halo mergers and is strictly related to 
the merging dynamic we want to estimate. Since we are ignor- 
ing this term, the resulting signal amplitude is smaller, making 
this choice conservative. In any case our filter can easily in- 
clude this component, providing an estimate of the dynamic of 
the whole infalling dark matter and not only of the main merger 
components. 



ILacev & Cold ill 993I) based on the extended Press-Schechter 
formalism to model the distribution of cluster mergers in the 
Universe. This description basically relies on the probability 
of a region of mass Mi having overdensity 61 provided that, 
at the same point, the overdensity measured on a larger scale 
(corresponding to a mass M2) is 62- In other words, our start- 
ing point is the conditional probability of having 61 on scales 
of Ml, given that on scales corresponding to M2 the density 
contrast is 62'- 



fi6uMi\62,M2)dSi 



■S2 



V2^(5i -82)^'^ 
(<5i - 62? 



exp 



2(5 1-52) 



(3) 



dS 1 



Here, S , refers to the variance of the 6 field on the scale of mass 
Mi, Si = {cr(Mi)^). Via the Bayes theorem, it is possible to 
compute the reverse probability, i.e., the probability of having 
a region of mass M2 with overdensity 62 given that on scale of 
Ml the density contrast is 6\. If one introduces now the rate of 
growth of density fluctuations in linear theory, then this reverse 



probability can be understood as the rate at which a given halo 
collects matter: 



'R{Mi,AM,z) = 



dhMdz 



nD{zfS2 



dD{z) 



dz 



do-{M2) 



exp 



dM2 

6l 



1 



(i-lf) 

1 1 
'S~2~'s'i 



3/2- 



(4) 



2D(z)2 

where, 6sc stands for the spherical collapse critical overdensity, 
and AM is the amount of mass being collected by the halo of 
mass Ml, so that M2 = Mi H- AM. The quantity D(z) is the 
linear growth factor of density perturbations. When compared 
to numerical simulations, this formalism is able to reproduce 
appro ximately the halo mer ger history (e.g. I\^ n den Boscll 
2002). However, as noted by Benson et alJ ( l2005l) . it has some 
intrinsic inconsistency due to its asymmetry in the arguments 
Ml and AM. In fact, the amount of clusters of mass Mi to which 
clusters of mass AM are merging does not equal the number of 
clusters of mass AM to which clusters of mass Mi are merging. 
We overcome this problem by defining the number of mergers 
of clusters of masses Mi and AM as: 



N^erg(Mu^M,z) = 



'R{MuAM,z)n(Muzh 



'R{AM,Muz)n(AM,z) 



(5) 



where n{M, z) re presents the nu mber density of haloes of mass 
M at redshift z, JPress & Schechteri.il9 74). Since Equation (0} 
provides a merger rate, we must multiply it by the finite 
redshift interval Az which is taken to be equal to the dynamical 
time of the final cluster (approxi mated by tdvniz) - Q.Q9/H(z) , 
with H(z) the Hubble constant ( lEnBlin & Rotteeringt l2002i) . 
Although this description of the halo merger rate might not be 
accurate, it should still provide a consistent description of the 
halo formation history. 



We follow the model given in, ISarazinI ( l2002t) . when as- 
signing distances and velocities to the merging clusters. The 
typical initial distance between a pair of clusters, following 
Equation (10) in Sarazin ( .200Z) . may be approximated as 



(6) 



where t{z) is the cosmic time. We assume that the probability of 
finding two merging clusters at a given distance r(Q < r < do) 
is directly proportional to the time spent by those clusters in the 
interval [r, r + dr], that is. 



P(r) <^ -7--, 
v(r) 



(7) 



where v(r) is the relative infall velocity of the clusters, given by 



v(r) ^ V() 



Ml H- AM 



Mo 




(8) 



In this equation, vq - 2930 km s ' and di -I Mpc, and we 
ignore the angular momentum of the merging system. Putting 
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all this together, we end up with a semi-analytical model of the 
merging cluster population in our universe, and of their veloci- 
ties. Finally, we assume an isotropic distribution for the orien- 
tation of the merger plane with respect to the observer. 

4. Observational strategy 

The millimetric observations of upcoming experiments will 
provide extensive catalogues of galaxy clusters detected 
through their thermal SZ effect (tSZ), together with estimates 
of their masses and their geometrical properties. The tSZ ef- 
fect can be measured separately from the other effects due to 
its unique spectral dependence. The data of CMB temperature 
fluctuations will contain signatures of the RS effect we aim 
to detect, together with the following components which will 
be regarded here as contaminants: the primordial CMB field, 
the instrumental noise, the lensing of the CMB anisotropics in- 
duced by the matter present along the line of sight, and the 
kinematic SZ (kSZ) eff'ect generated by the cluster peculiar ve- 
locity with respect to the Hubble flow. As long as the contami- 
nants show different spatial characteristics compared to the sig- 
nal we are trying to unveil (the RS effect), it is possible to define 
a filter which optimally reduces the impact of those sources of 
noise. Regarding this, instrumental noise will be assumed to 
be scale-free white noise, i.e., to have the same power at all 
scales. The intrinsic CMB fluctuations usually have (if clusters 
are not too close) scales larger than those of the RS signal. On 
the other hand, the kSZ and the tSZ residuals are confined to 
the region of the cluster where the gas is located. Compared to 
these scales, the RS is usually broader 

However, if a contaminant shows spatial power on scales 
which are too close to those of the signal, only a more direct 
approach can be performed. This turns out to be the case for 
the anisotropics introduced by the lensing of the intrinsic CMB 
fluctuations. This component will require a separate processing 
(Section lOl : from a template of the lensing-induced deflection 
field, we shall "deflect back" the data, in an attempt to subtract 
the effect of gravitational lensing on the CMB. This will allow 
us to minimize the impact of this contaminant. 

Once the cluster lensing signal is reduced by this "de- 
lensing procedure", an optimized filter will be applied to sup- 
press the other noise terms and to estimate the RS effect am- 
plitude. This amplitude depends on the cluster infall velocity, 
which is essentially the physical parameter we are measuring 
(see Equation |2}. This filter incorporates in its definition two 
templates: one for the RS effect described in Section pa- 
rameterized by the mass which we assumed to be estimated 
from observations of the tSZ effect, and one for the kSZ effect. 
The second one could be directly derived from the millimet- 
ric observations by observing at 217GHz and filtering away the 
CMB features through a high-pass filter 

As we shall see, the resulting signal-to-noise ratio of a sin- 
gle merger is too small to provide a detectable signal. Thus it is 
necessary to average multiple measurements of many mergers. 
The detection of the cluster throug h their tSZ effec t is akeady 
well described in other papers (e.g. lSchaefer et al 112004) so It 
will not be discussed here. We will now describe in detail the 
steps of the observational strategy summarized above. 



4.1. De-lensing procedure 

The CMB anisotropy produced by the gravitational lensing ef- 
fect of galaxy clusters has a typical amplitude ten times larger 
than the RS anisotropy induced by the same clusters. These 
two effects arise from the same gravitational potential, so that 
their typical spectral scales are comparable and correlated in 
position and, although their patterns are not identical'*, the sim- 
ilarity of their power spectra complicates the use of any Wiener 
filter 

Therefore it is necessary to reduce the lensing effect before 
the application of any optimized filter To achieve this, we pro- 
pose to remap the observed data by applying a distortion which 
compensates, as accurately as possible, the cluster lensing ef- 
fect. 

To do this it is necessary to invert the lens equation T(0) - 
f(6 - a(6)) by introducing the new coordinate 0' - 6 — a{6). 
In this way it is possible to derive from the observed data T 
the de-lensed map f thanks to the deflection field template a, 
according to 

f (61) = T(0 + a(0')) . (9) 

Note that the deflection field has to be computed via an iterative 
procedure, because it is evaluated at the new coordinate 6'. 

In order to build a deflection field template, we define the 
merger using two NFW haloes according to Equation Q. We 
assume that their masses are derived by the tSZ measurements, 
and for this reason we will consider some uncertainty in the 
mass determination. The assumption of spherical haloes is jus- 
tified since the deflection field is roughly symmetric on large 
scales. Furthemore, the final estimate of the signal will be pro- 
duced after stacking all available merger events. 

Of course this procedure will - apart from the CMB - 
also distort the RS signal, the kSZ signal and the distribution 
of the instrumental noise. Because of this, the RS effect and 
the kSZ templates used in the optimal filter definition have to 
be corrected through the same process. This procedure leaves 
some randomly distributed residuals in the field, and some non- 
Gaussian features in the instrumental noise which will be ig- 
nored in the following analysis. 

4.2. Filter 

Once the cluster lensing contribution is minimized, the re- 
maining data contaminants in the maps are the CMB primary 
anisotropy, the kSZ effect and the instrumental noise. 

The Gaussian contribution of CMB and instrumental noise 
are fully described by their corresponding power spectra and 
can be filtered out by optimized filters. In particular, the CMB 
power spectrum shows an exponential cutoff around the ar- 
cminute scale due to Silk damping of primordial fluctuations. 
This angular scale is close to the typical one for distant galaxy 
clusters, and hence the CMB intrinsic fluctuations should make 

^ The difference in their patterns arises from the fact that the RS 
effect orientation depends on the velocity direction of clusters, while 
the gravitational lensing anisotropy depends on the local CMB tem- 
perature gradient direction (see Section|2J. 
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Table 1. Main characteristics of the upcoming millimetric ob- 
servatories ACT, SPT and of our simulated instrument. The 
second column gives the available frequency channels, the third 
column the angular resolution and the last column the relative 
sensitivity. We are reporting nominal sensitivities for ACT, and 
those of SPT correspond to one hour of observation. 





bands 


FWHM 


AT/beam 




(GHz) 


o 




ACT 


145, 225, 265 


1.7, 1.1, 0.93 


2, 3.3, 4.7 


SPT 


150, 219, 274 


1,0.69, 0.56 


2.5, 3.0, 2.65 


simulation 


217 


1 


1 



small contributions at the spatial frequencies we are interested 
in. For CMB experiments with very high-sensitivity on small- 
scale, such as ACT, SPT or the Atacama Large Millimeter 
Array^ we shall assume that the noise is spatially white, i.e., 
that it introduces the same amount of power in all scale. 

We present here an optimally matched filter (derived de- 
tail in AppendixlAl. which maximizes the signal-to-noise ratio 
by processing the signal in the F ourier domain. A similar filter 
construction was proposed by .Haehnelt & Tegmark. (1996i) in 
order to extract the kSZ signal from single clusters. However, 
that method considered the CMB and instrumental noise the 
only sources of confusion when trying to recover the kSZ sig- 
nal. Our filter is defined in a slightly more general context. It 
provides an estimate of the amplitude of a signal with a known 
spatial template and is to be measured despite a presence of ho- 
mogeneous noise with a known spatial power spectrum and an 
arbitrary number of other contaminants following known spa- 
tial patterns. Hence, the total observed data can be modeled by 



s{r) = AT(r) + n(r) + ^ vmir) , 



(10) 



where A is the amplitude of the signal, T(r) is the signal shape, 
n(r) is a Gaussian random noise with a known power spec- 
trum Pn(k), and z,(r) are M noise sources with an amplitude 
V,. Note that this model can a priori be applied to different 
problems. For instance, in weak lensing observations, where 
it is necessary to take into account the noise introduced by 
stars, bright galaxies and the cluster core itself, this filter al- 
lows improved results obtained with standard matched filters 
(e.g. Maturi et al., 2005b). 

The Gaussian noise contribution «(r) will contain the CMB 
temperature fluctuations, the instrumental noise, and all the 
noise sources which can be modeled as Gaussian random fields. 
In this case, the filter derived in Appendix (|AJ takes the form 



AT(k) - BiZi(k) - B2Z2ik) 
P„(k) 



(11) 



where Pn(k) - PcMsik) + Pinsi.(k) is the total Gaussian noise 
power spectrum, including the CMB and the instrumental 
noise. In our case, we introduce two spatially characterised 
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Fig. 2. To give an intuitive view of how our filter works, we 
present here two simple examples where the signal and a con- 
taminant of known spatial pattern have both Gaussian profiles 
with the same width. In the upper panel both signal and noise 
are positive and the sum of these contributions, the data, is 
given by the solid line, while the dashed line displays the re- 
sulting filter In the bottom panel we only changed the sign of 
the amplitude of the noise, which now is negative. The esti- 
mator defined in Equation dlSI l is the integral of the solid line 
times the dashed one. In the first (upper panel) case, the filter is 
negative where the contaminant is positive, so that their prod- 
uct is negative and the contaminant contribution is subtracted. 
In the second case (bottom panel), both filter and contaminant 
are negative so that their product is positive and the negative 
contribution of the noise is compensated. This actually shows 
that the same filter provides an unbiased estimator also if only 
the variance of the contaminant amplitude is known. That is, 
the sign of the contaminant is irrelevant to the filter, and this 
reflects the statistical nature of our filter. 



contaminants (M - 2), that is, the kSZ signal of the two merg- 
ing haloes Z[(k) and Ziik). The normalization of the first term 
in the numerator can be expressed as /I = (1 + D)/C, with 



ALMA, see www.eso.org/projects/alma 



C = [\Tf/Pn} 

and 

^_^^, r(BiZ,+B2Z2) \ 



The coeflicients B, are given by 
Bi^crl{Zy}, 



(12) 



(13) 



(14) 
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Fig. 3. Spatial distribution of the signals considered in a ma- 
jor merger at z = 0.3, of two haloes with masses lO'^M© and 
5 X IO'^^Mq, and with an infall velocity of 1500km/ s at 45 de- 
grees with respect to the line of sight. Since the two haloes 
are colliding, the RS effect (top left) of both components adds 
constructively, amplifying the resulting signal in the region be- 
tween them. The kSZ effect (top right) shows a typical dipo- 
lar pattern reflecting opposite velocities along the line of sight. 
The anisotropy introduced by the lensing effect has a complex 
shape, since it depends on the configuration of both haloes and 
the CMB gradient (bottom left). The bottom right panel shows 
our filter It is high-pass in order to suppress the CMB primary 
anisotropy noise contribution. The temperature amplitudes of 
these panels along the line containing the haloes are displayed 
in Figure 0. 



and the estimate of the RS effect amplitude A will be provided 
by 



A = {S'¥*] 



(15) 



In these equations, the bracket notation { / ) stands for a double 
integral in the k-space of a general function f{k) as defined 
in Equation (IA.2> . The quantity cr^. represents the variance of 
the ith contaminant amplitude, and the symbol % denotes "real 
part". 

It is clear from Equation (II H that the filter can be spUt into 
two components. The first, (T* /P„), maximizes the sensitivity 
at the spatial frequencies where the signal T is large and the 
noise power spectrum P„ is small, as opposed to the second 
{-(B\Zi + B2Z2)/P„), which underweights the regions where 
the Z, components are important and correlated with the filter 
itself. Note that since the filter is present in the definition of 
the Bi coefficients, the final form of the filter must be computed 
iteratively. 



The template for the RS effect (T(r)) was modeled by two 
merging NEW haloes, as explained in Section (|2jl. The con- 
tribution from small merging substructures should be negligi- 
ble, since our filter is a weighted integral over the whole field 
(Equation 1 15l l. and these minor mergers are odd functions of 
typical scales much smaller than the one of the filter 

The kSZ masks are applied as a weight templates and are 
not directly subtracted from the input data. This makes the ac- 
curacy to which the templates are known uncritical. Eor this 
reason it is possible to use the observations at 217 GHz (the 
frequency at which the tSZ vanishes to a good approximation) 
as a good template for the kSZ. There might be some tSZ resid- 
uals, together with the lensing and RS signals, but their contri- 
bution is negligible for the filter construction. Again, this is 
because the filter does not subtract this template directly from 
the maps, which of course would erase all the RS signal we are 
searching for, but it only weights down optimally the regions 
contaminated by the kSZ effect. 

5. Simulation of the observations 

To test the described procedure we performed simulations of 
the secondary anisotropics induced on the CMB by binary clus- 
ter mergers in a cosmological context. Eor the two merging 
haloes, we adopt a NEW profile for the dark matter compo- 
nent and a spherical isothermal sphere (SIS) for the gas distri- 
bution. We ignore tidal effects and the hydrodynamic processes 
involved in the interaction between the intra-cluster gas (ICM) 
of the two haloes. Eor our purpose, the model adopted gives a 
sufficiently reaUstic representation of RS, kSZ, lensing effect 
and instrumental noise. 

We adopt a standard ACDM cosmology with a density con- 
tribution from dark matter, baryons and cosmological constant 
of Qdm = 0.27, Q.B = 0.044, and Oa = 0.7, respectively. 
The Hu bble constant was set to Hq = 100 h km s"' Mpc"' with 
h = 0.7 jBennett et alll2003l) . 

We simulated the CMB primary anisotropics as Gaussian 
random fields with a resolution of 1024^ pixels and a field of 
view of 120'. The CMB power sp ectrum used was computed 
using CMBEASY ("Ooraiil l2003l) . assuming a re-ionization 
fraction of 0.1 at redshift z = 6.2 and a present helium abun- 
dance of 0.24. The field size ensures an adequate sampling of 
the CMB multipoles, and permits the retention of only the cen- 
tral quarter of the whole field, avoiding the boundary effects of 
the Fourier transforms. 

For an isothermal cluster with negligible internal motions 
the kSZ and the tSZ effects are proportional to one another: 



-T^W = 8tsz(y) — 2 



dl rieCTT, 



Tq 



(0) 



dl rieCTT. 



(16) 



in the previous equation gtsziv) provides the (non-relativistic) 
frequency dependence of the tSZ effect in temperature, and nig. 
He and Te are the electron mass, number density and tempera- 
ture, respectively. The integrals are performed along the line of 
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Fig. 4. Comparison between tfie filter profile and the different components under consideration. The top left panel shows the 
de-lensed RS effect superimposed on the filter. Note the sudden change of amplitude and sign of the filter close to the cluster 
cores. The top right panel displays the kSZ effect and the filter, whereas the bottom left one shows how the de-lensing procedure 
acts on the kSZ effect, shrinking its profile and shifting the peaks closer to the center In the bottom left panel we show how 
the de-lensing procedure acts on the kSZ effect, by shrinking its profile and sh ifting the peaks closer to the center. The bottom 
right panel compares our filter (new filter) with the filter originally proposed bv LSeliak & Zaldarriagal l EoQol) (H & T 1996). The 
difference between these two filters is negUgible when the contaminant given by kSZ effect (the Z template) is orthogonal to the 
signal, otherwise it becomes relevant. 



sight given by the pointing direction n, and the angle 6 sepa- 
rates the line of sight from the cluster center. The dependence 
of the integral on the line of sight translates into the angular 
dependence. In particular, we assumed the ICM of each merger 
component to be distributed as an isothermal sphere with a ra- 
dial y6-profile with /3 - I (i.e. a King profile). We based our 
strategy on observations at 217 GHz, where the tSZ effect van- 
ishes. Thus the SZ contamination is given only by the kSZ com- 
ponent. 

The lensing effect on the CMB is described by the lens 
equation: 



T(0) = f [61 - a(e)] ~ T(0) - VT(e) ■ a(0) , 



(17) 



where T is the observed temperature map, t is the un - lensed 
CMB and a is the cluster deflection angle. iBartehnannl i 1 996l) 
wrote an analytic expression for the deflection angle caused by 
a NFW halo, and, by using it we can compute the VJ ■ a term. 



which rep resents the secondary anisotr opy given by the lens- 
ing effect (Se hak & Zald arriaga. '200(f). The linear expansion 
is justified as long as the variations of the local CMB tempera- 
ture gradient are small on a deflection angle scale. 

The RS and kSZ effects were computed after introducing 
this deflection angle into Equation 0. The distribution of clus- 
ter mergers in the universe is described according to the equa- 
tions presented in Section 0. The number of clusters Nm,z 
was computed with that formalism on a finite number of mass- 
redshift cells. Then, for each cell, we randomly distributed Nm,z 
clusters with redshifts and masses enclosed in the cell range. 

We included in our simulation the noise from a single dish 
telescope with a resolution of 1 arcminute and a sensitivity of 
1 fiK per beam. This noise was modeled as a Gaussian field with 
the following power spectrum 



Cf"^'= = exp 



1(1+ 1)FWHM2 



81n2 



(18) 
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where AT IT is the sensitivity of the experim ent in on the 
scale of the beam and w ' (Ar/rFWHM)^ ( lKnoxL [l995). 
The exponential expresses the increasing weight of noise when 
probing sub-beam scales. These characteristics of the simulated 
instrument are compatible with the next-generation millimetric 
observatories, as shown in Table Q. We have not included the 
parameters corresponding to ALMA because our filter should 
have then included the synthesized beam of the interferometer, 
but its sensitivity should enable it to detect signals of ampli- 
tudes similar and smaller than those of the RS effect considered 
here. 



6. Results 

In order to probe the ability of our extended filter in separat- 
ing the RS signal, we have prepared three different mock cat- 
alogues of cluster mergers following the recipe described in 
Section @. These catalogues were generated using three mini- 
mum mass thresholds of 5 x lO'^Mo, 8 x IO'^Mq and W^M^, 
and contain 6807, 266 and 28 cluster mergers, respectively. 

The analysis of the mock catalogues was carried out in two 
ways, corresponding to the degree of knowledge of the tem- 
plate maps. 

- Case A. We assume that we know the masses of the merg- 
ing clusters with an uncertainty of 30%. These values are 
used to build the kSZ and the RS templates. For the halo 
radial velocity, we use Equation MA\ with an expected ve- 
locity dispersion cr,. k 300km/ s. 

- Case B. Similar to the previous case but, instead of build- 
ing a kSZ template, we shall use directly the millimet- 
ric observations of the kSZ effect, which ac count for the 
kSZ asymmetries caused by internal motions ( Nagai et'all 
I2OO3.) . The presence of the RS effect in these maps is neg- 
ligible and does not affect the filter definition. The only 
problem related to this procedure is the noise present in 
the kSZ maps which is included in the filter construction. 
A low pass filter of the map should reduce the problem. In 
our simulation we directly used the kSZ plus the RS maps, 
assuming that a noise suppression procedure had akeady 
been applied to those maps. 

The results for Case A are presented in Table Q. The 
columns represent: first, the lower limit adopted to build the 
cluster samples; second, the number of objects in the sample; 
third, the average RS signal of the sample; and last, the esti- 
mate. 

Even with a moderate size cluster merger catalogue (266 
clusters), the significance of our method approaches the 2-cr 
level in the detection of the RS signal. Since we are observing 
different mergers, and the errors in different mergers are un- 
correlated, our Poissonian error bars scale is the inverse square 
root of the number of mergers. 

The results for Case B are also presented in Table As 
one can see, RS signal detection is still possible. Comparing 
these results with those of the previous table, it is clear that 
the filter efficiency is not significantly affected. This is mainly 
due to the statistical nature of our filter: it rehes not on the 



Table 2. RS amplitude estimates for three mass selected cluster 
subsamples, defining the filter with a 30% mass uncertainty and 
using only the expected velocity dispersion cr,, x 300 km/s. 
Error bars in the last two columns correspond to the combined 
measurement from all clusters in the corresponding sample and 
with the coresponding filters (case A and B). 



mass 


clusters 


simulation 


estimate A 


estimate B 


Me 




fiK 






M > 5 X 10'-' 


6807 


0.18 


0.13 ±0.03 


0.16 ±0.03 


M > 8 X 10''' 


266 


0.25 


0.24 ± 0.14 


0.13 ±0.19 


M > 1 X 10"* 


28 


0.28 


0.29 ± 0.42 


0.05 ± 0.68 



accuracy of our templates for each merger, but on the accuracy 
of our templates when describing the average spatial properties 
of the kSZ and RS signals. Indeed, in principle one could also 
use the same average RS and kSZ templates for all mergers. 
The precision of our method is hence limited by the number of 
available mergers according to Poisson statistics. 

We recall that these results have been obtained ignoring the 
component of the RS signal of galaxy clusters given by the 
collapse of the environmental dark-matter in the potential well 
of the system. The inclusion of this component would increase 
the signal-to-noise ratio of the estimates. 

7. Conclusions 

We presented a method to extract the kinematic properties 
of major merger events of galaxy clusters by means of their 
RS signal, which is observable at centimeter/millimeter wave- 
lengths. The RS effect is a secondary anisotropy of the CMB 
produced by the time variation of the gravitational potential 
along the line of sight. In this particular case, the main contribu- 
tion comes from the infall motion of the merging components, 
so that the total signal is a measurement of the convergence 
of the projected perpendicular momentum of the system (see 
Paper I). 

We proposed a method to extract this RS signal from a 
sample of cluster mergers observed at millimetric wavelengths. 
These observations will be contaminated by the primordial 
CMB fluctuations, the kinematic SZ signal and the lensing ef- 
fect from the clusters, together with the instrumental noise. 
We assumed that the thermal component of the SZ effect can 
be subtracted by multi-frequency observations, given that we 
know its frequency dependence. 

The RS effect amplitude in one single merger event is much 
smaller than all the other components and thus only a statisti- 
cal detection of the signal may be achieved. To maximize the 
signal-to-noise ratio and decrease the number of mergers that 
we need to co-add, we can apply a filter which enhances the RS 
signal above the other components. Therefore, we proposed in 
Section (I4.2jl an extended versio n of a matched filter presented 
by Haehn elt & TegmarkI ( Il996l) . A complete derivation of our 
filter is given in Appendix (|A|i. It requires as input the power 
spectra of the CMB and the instrumental noise, plus two tem- 
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plates, one for the RS signal and another for the contaminat- 
ing kSZ effect. These templates can be simple analytic mod- 
els parameterized with the mass estimates provided by the tSZ 
measurements, and with the expected average velocities. We 
would like to note that the filter construction is general so that it 
can be applied to other problems l ike the extraction of gravita- 
tional lensing on CMB maps (e.g. lHaehnelt & T egmark. 1996|: 
IMaturi et al., 2005a|) or on background galaxies (Maturi et all 

We showed that the developed filter can control the contam- 
ination from all signals except CMB lensing. This is because 
the RS and the lensing effect originate from the same gravi- 
tational deflection field, and thus their spatial frequencies are 
similar Therefore, we included in our pipeline an additional 
step before the filter application, i.e. a "de-lensing" procedure. 
We propose remapping the input data according to the inverse 
template deflection field, again parameterized by the tSZ mass 
estimate. 

The method was applied to mock catalogues of cluster 
mergers. Our results show that the RS signal of merging clus- 
ters could be measured with this method after observing of 
the order of 1,000 cluster mergers. These numbers will be 
achieved by combining the expected yields of the upcoming 
high-resolution millimetric surveys (e.g. ACT, SPT, ALMA). 
Assuming a fraction of 30% for the number of mergers over the 
number of clusters, a RS detection will be feasible if future SZ 
surveys are able to perform sensitive observations of the order 
of 10"* clusters. 



Appendix A: Optimal filter 

We aim to build an optimized filter iff{r) to extract from a noisy 
data set s{r) the best estimate of amplitude A of any signal with 
a known spatial shape T(r): 



s(r) - AT{r) + n{r) + z(r) . 



(A.l) 



We assume n(r) to be a Gaussian random noise with zero mean 
({n(r)) - 0) and known power spectrum PNik); and z(r) = 
Yjfii ViZi(r) to be the sum of M different noise components with 
unknown amplitude v, with zero mean ((v,) - 0) but known 
variance cr,,,-. 

We define the integral of a general function f{k) as 



(A.2) 



Since the Equation (IA.1> is linear in A, its most general 
linear estimate in the Fourier domain can be written as 



(A.3) 



where the capital letters refer to Fourier transforms. This esti- 
mate is required to have a bias b - {A -A) - 0. It has a variance 
cr^ - {{A - A)^) given by 



/7=A({r>F*)- 1) 



(A.4) 



where we grouped the collection of the M noise contributions 
Zi{k) in a vector Z{k) and defined the covariance matrix V of 
its components as Vij - (v,v'j). 

Since we are interested deriving a filter which minimizes 
the estimate variance maintaining the unbias condition b - Q, 
we introduce the Lagrangian multiplier A and search for the 
filter function which minimizes the action L - cP- + Ab, ob- 
taining 



AT(k)- Z'(k)V{Z*'i'} 



(A.6) 



where the apex f is for the transposed vectors. 

The normalization factor A is obtained by substituting 
Equation iA.6\ in Equation JA.4> . yielding 



A 



1+%([T*Z'/P„}V{Z*^'}) 

myPn] 



(A.7) 



where is for the real part. 

As is clear from Equation iA.6\ . the filter has to be com- 
puted through an iterative procedure, where the starting point 
can be the filter with the variance matrix V equal to the null 
matrix (this is essentially the filter originally proposed by 
'Haeh nelt & TegmarkI \l996)). This iterative procedure is re- 
quired since the filter evaluation requires the correlation be- 
tween the filter itself and the noise sources. 

The simple intuitive interpretation of this filter is that the 
first term T/P,, maximizes the sensitivity on the spatial fre- 
quencies where the signal T is large and the noise power 
spectrum P„ is small. At the same time the second term 
{Z'V{Z*^'])/P„ introduces a weighted mask which takes into 
account the correlation between the filter *P and all the Z, noise 
components, as well as the correlation between the different Z, 
noise components. It is important to note that the filter 'I' does 
not apply a mask only by subtracting the Z, templates from the 
input data, leading to spurious residuals, but it down-weights 
the influence of the regions where the Z, components are dom- 
inant and correlated with the filter itself. This means it is not 
critical to have a very accurate template for the noise sources, 
in order to apply the filter properly. It is sufficient to have a 
simple model which represents their average shapes. 

The derived filter satisfies a generalized problem, where the 
noise sources v,Z, can also be correlated, but in the majority of 
cases the noise sources will be uncorrected, either due to them 
having different positions or their amplitudes v being mutually 
uncorrected. In this case the matrix V is diagonal and the filter 
function simplifies to 



AT(k)-Y.fBiZi(k) 
Pnik) 



(A.8) 



where the normalization can be written as /I - (1 + D)/C, the 
denominator C is given by 



C = {\Tf/P„) 

and the Z?,'s and D are defined as 



(A.5) Bi^irUz;^']. 



(A.9) 



(A.10) 
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and 

^^^^^ 



With this assumption, valid for many applications, the filter ac- 
quires this simple expression. 
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